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ABSTRACT 


Heat-transfer, thermal bending, and mechanical buckling analyses have been performed on a superalloy 
“honeycomb” thermal protection system (TPS) for future hypersonic flight vehicles. The studies focus on the 
effect of honeycomb cell geometry on the TPS heat-shielding performance, honeycomb cell wall buckling 
characteristics, and the effect of boundary conditions on the TPS thermal bending behavior. 

The results of the study show that the heat- shielding performance of a TPS panel is very sensitive to change in 
honeycomb core depth, but insensitive to change in honeycomb cell cross-sectional shape. The thermal 
deformations and thermal stresses in the TPS panel are found to be very sensitive to the edge support conditions. 
Slight corrugation of the honeycomb cell walls can greatly increase their buckling strength. 

NOMENCLATURE 

a depth of honeycomb core, or length of a rectangular plate, in. 

a' = la, 1.5a, 2a, modified depth of honeycomb core, in. 

(0 9 

A ( , cross-sectional area (perpendicular to the cell axis) of honeycomb cell type i, in“ 

(0 9 

A w cross-sectional area (perpendicular to the cell axis) of cell walls of honeycomb cell type i, in 

b width of a rectangular plate, in. 

bj length of free side of type i honeycomb cell cross section, in. 

b' length of bonded double side of TPS cell cross section, in. 

b = 5 c, corrugation length of corrugated side of corrugated cell cross section, in. 

c length of one corrugation leg, in. 

dj size of type i honeycomb cell (maximum diagonal of the cell cross section), in. 

Dj bending stiffness of wall panel of type i honeycomb cell, in-lb 

E Young’s modulus, lb/in“ 

Fjj radiation view factor, defined as the fraction of radiant heat leaving radiation element i and incident on 

radiation element j 

G shear modulus, lb/in~ 

h c corrugation depth (twice the amplitude of corrugation wave), in. 

i index (1, 2, 3,...) associated with type i honeycomb cell geometry: i = 1 for right hexagonal cell; / = 2 

for square cell; i = 3 for TPS flat cell walls; i = 4 for TPS corrugated cell walls 



i index (1, 2, 3,...) associated with radiation element i 

j index (1, 2, 3,...) associated with radiation element j 

JLOC joint location 

,2 

b- n\ 

k = — — (N z ) c ' , buckling coefficient for a rectangular plate with aspect ratio a/b- under uniaxial 

ft D i compression in a direction 

/ side length of square honeycomb-core sandwich TPS panel, in. 

m number of buckled half-waves in rectangular plate longitudinal direction 

n number of buckled half-waves in rectangular plate lateral direction 

(/V classical compressive buckling load intensity in the z direction of the weakest rectangular wall panel of 

type i honeycomb cell, lb/in 

l~2 2 

p one-half of corrugation pitch, Jc~ - h ( . , in. 

q aerodynamic pressure over outer surface of honeycomb-core sandwich TPS panel, lb/in“ 

( f) ^ 

q cr critical value of q at buckling for type i honeycomb cell walls, lb/in“ 

SPAR Structural Performance and Resizing (finite-element computer program) 

STS space transportation system 

t time, sec 

t honeycomb cell wall thickness, in. 

~t coiTugated honeycomb cell wall thickness, in. 

t sandwich face sheet thickness, in. 

T , inner face sheet temperature, °F 

TPS thermal protection system 

T(t) transient surface temperature over outer surface of honeycomb-core sandwich TPS panel, °F 

T outer face sheet temperature, °F 

v Poisson’s ratio 

w g deflection at sandwich panel edge midpoint, in. 
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w 

max 

w 


x, y, z 


a 


AT 

6 

Pc 


Pi 



CF 


cr 


cr_ 


(<r/> 

z' cr 


Subscripts: 


maximum deflection at sandwich panel center, in. 
weight, lb 

rectangular Cartesian coordinates 
coefficient of thermal expansion, in/in-°F 

= T - T i , difference between outer and inner face sheet temperatures, °F 
corrugation angle, deg 

density of material used for honeycomb core and sandwich face sheets, lb/in 

o 

effective density of honeycomb core with type i cell geometry, lb/in 
effective density of sandwich panel with type i honeycomb cell geometry, lb/in 
stress in x direction induced in the face sheets, lb/in“ 
stress in y direction induced in the face sheets, lb/in“ 

compressive stress in z direction in the honeycomb cell walls induced by the surface pressure q 
compressive stress in z direction in the honeycomb cell walls at the buckling of the weakest cell wall 


( ) critical value at buckling 

( ) • associated with type i honeycomb cell geometry 

Superscripts: 

( )^ associated with type i honeycomb cell geometry 


INTRODUCTION 

Flypersonic flight vehicles such as the Space Shuttle orbiter are subjected to severe aerodynamic heating 
during flight missions. A thermal protection system (TPS) made of low-thermal conductivity materials is used to 
insulate primary structures from overheating so that the vehicle can operate within the design temperature limit 
(ref. 1). Typical types of TPS include porous silica rigid “tiles” (ref. 2), flexible felt “blanket” (ref. 3), multiwalled 
titanium panels with dimpled cores (refs. 4, 5), and so forth. The main disadvantage of the rigid TPS tiles is the 
fragility of the materials used. After each mission, numerous tiles must be replaced because of surface damage 
incurred. To eliminate such repairs, a superalloy “honeycomb” TPS is currently being considered for future space 
transportation systems (STSs) (refs. 6, 7). The honeycomb construction provides low density and low thermal 
conductivity through the TPS thickness. The superalloy TPS is capable of functioning at high temperatures (ref. 8). 
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A TPS is exposed to high temperatures on the outer surface and to relatively lower temperatures on the inner 
surface facing the cooler substructures. This thermal gradient across the TPS thickness will induce thermal bending 
moment, causing the TPS panel to bulge outwardly. A highly efficient TPS could provide large temperature 
gradient across its depth; however, excess temperature gradient across the depth could induce excess thermal 
bending of the TPS panel, thus disturbing the local air flow. Therefore, proper trade studies must be exercised to 
design optimum honeycomb TPS structures. Also, the honeycomb TPS will be subjected to aerodynamic pressure 
loading, and cell walls must be stiff enough to withstand the air load without buckling. 

This report presents the results of the heat- shielding and thermostructural analyses of honeycomb “sandwich” 
TPS panels that have different honeycomb cell geometries. 


DESCRIPTION OF PROBLEM 


Figure 1 shows a honeycomb-core sandwich TPS panel subjected to transient surface temperature, T(t ) , and 
aerodynamic pressure, q, over its entire outer surface. The TPS panel is square with a side length of /, and is 
fabricated with two identical face sheets with a thickness of t and a honeycomb core with a depth of a. For a 
given material, the overall heat-shielding and thermostructural performances of the honeycomb TPS panel depend 
on the thickness of the face sheets, depth of the honeycomb core, thickness of the honeycomb cell walls, and size 
and shape of the honeycomb cells. 

The first problem is to perform geometrical (shape) analysis of different candidate honeycomb cells that have 
the same effective density but different geometrical shapes. The second problem is to perform heat-transfer 
analysis of TPSs with different honeycomb cell geometry. The third problem is to perform the thermal bending 
analysis of the TPS panels subjected to one-sided heating and under different support conditions. The fourth 
problem is to perform comparative buckling analysis of different types of honeycomb cell wall panels under 
compressive surface air pressure. 


HONEYCOMB STRUCTURES 


This section presents the geometrical analysis of honeycomb cells with different geometry. The honeycomb 
TPS fabrication process also is described. 


Honeycomb Construction 

The honeycomb-core structure is fabricated using multiple layers of thin strips of metallic (or nonmetallic) 
foils joined together and properly deformed (ref. 9). The thin strips are bonded together first at equally spaced 
parallel zones (fig. 2). The equally spaced bonding zones on one side of each thin strip are staggered with respect to 
those on the opposite side. The bonded multiple sheet assembly then is pulled apart in the thickness direction 
through bending of the bonded and free junctures to form a final honeycomb structure (fig. 2). By modifying the 
width of the interfacial bonding belt zones and their spacing, and by pulling apart the multiple layers through 
bending deformations to a desired level, a family of different honeycomb cell geometry could be generated. 
Figure 3 shows the exploded view of the construction of a sandwich panel with upper and lower face sheets to be 
bonded to the honeycomb core through the brazing process. 
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Honeycomb Cells 


Figure 4 shows four types of honeycomb cell geometry to be analyzed. The honeycomb cell wall thickness for 
the first three types is t . The first type is a right hexagonal cell (fig. 4(a)) with identical side lengths of b^ . The 
second type is a square cell (fig. 4(b)) with side lengths of b 2 , which is modified from the right hexagonal cell by 
reducing the bonding interface length to a minimum of J2 1 . The third type is a distorted honeycomb cell with flat 
cell walls (fig. 4(c)). In this type of honeycomb cell, the bonded side length b' is slightly less than the free side 
length Z? 3 , and the angle between the two adjacent free cell walls is adjusted to be 90 degrees. 

The fourth geometry type is also a distorted honeycomb cell, but the free cell walls are corrugated (fig. 4(d)). 
Each corrugated wall has five corrugation legs. Cell wall thickness for the fourth type is either t ( , or t c (</.. ), 
which decreases with increasing corrugation depth, h c , for maintaining constant core density (fig. 5). The third and 
fourth types of distorted honeycomb cells have the same linear dimensions (that is, identical b' and Z? 4 = Z? 3 ), and 
are the approximate cell geometry of the proposed superalloy TPS for a new STS. The size, d-(i = 1, 2, 3, 4) , of 
each type of honeycomb cell is defined as the maximum diagonal of the cell cross section (fig. 4). 

Effective Densities 

If and , respectively, are the cross-sectional areas (perpendicular to the cell axis) of honeycomb cell 
type i and its cell walls, then the effective density of honeycomb core with type i cell geometry, p,- (face sheets 
excluded), is given by the following: 


p 

- = ; (/= 1,2, 3, 4) 

A l > 


(l) 


where p c is the density of the material used for the honeycomb core and sandwich face sheets. 

The density of the sandwich panel (including honeycomb core and face sheets) can be written as the following: 


o (0 A {t) a + 2 1 A (i) 

Psw _ a + n s A c 

Pc A {l) a + 2tA { ! ] 


Pc a _ 

2t s 

1 + — 

a 


( 2 ) 


in which equation (1) was used. The weight W of the TPS panel is given by 


W - p { 2f(a + 20 = Vp ( 


\ 


it ^ P ‘ 

2 1 + a — 

Pc, 


( 3 ) 
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In the calculations of , only one-half thickness of the cell wall was used because each of the cell wall 
panels is mutually shared by two adjacent cells. Namely, the unit cell wall thickness is 2t ( ./2 = t for the bonded 
double walls and t / 2 for the free walls. The effective densities for the four types of honeycomb cores can be 
written as follows (see the appendix for detailed calculations): 

Right hexagonal cell (i - 1): 


Square cell (i = 2): 


Pl _ _ Zj3t c _ I6.j3, c _ _ d, 

Pc A [ l) 9b i 9d l ' ' 2 



c 


TPS cell (flat wall) (/ = 3): 


P3 

Pc 



1(3 + V2)p 

5 b 3 


8^(3 

15 v 




( 4 ) 


( 5 ) 


(6) 


TPS cell (coiTugated wall) (/ = 4): 


Pa 

Pc 


OJ ' 
1 53-1 

+ 

oo 

II 

( 

3 

1+ 

1 b A ) 

b 4 15 

i *1 

3 4 J 



in which the following expression for coiTugation length b was used (fig. 5 and appendix): 


b_ 

b A 


1 + 24— 

\ b l 


l + 2J6K 

3 4 


( 7 ) 


( 8 ) 


Cell Sizes 


To compare the heat-shielding and buckling performances, the size of honeycomb cells types 1, 2, and 4 are 
adjusted to have the same effective density as that of honeycomb cell type 3 (that is, p 1 = p 1 = p 3 = p 4 ). By 
equating equations (4) and (6), and equations (5) and (6), the cell sizes d j and d 0 are determined as 


d i 


10 ^ d, 

3 ^ 2(72 + 3 ) 


( 9 ) 
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15 


(10) 


d 2 _ 


4(72 + 3 ) 


Figure 6 shows the plots of d 1 and d 2 as functions of d 2 based on equations (9) and (10). For the corrugated 
cell, by definition 


d ^ — d 3 


( 11 ) 


and to maintain the same effective density (that is, p 4 = p 3 ), the corrugated cell wall thickness t ( , is reduced as 
the corrugation depth h ( . increases according to the following equation (see appendix). 


t 

t 


c 

c 


3 + 72 


3 


1 + 24— + 

^ b l 


72 


3 + 72 


|i + ^ + 72 

3 4 



(12) 


Figure 7 shows the corrugation angle, 0, plotted as a function of corrugation depth h/t ; figure 8 shows the 
corrugated cell wall thickness t c /t decreases with the increasing corrugation depth h /t . 


Numerical Input Values 

A typical candidate TPS panel for a new STS has the following dimensions: 

a = 0.488 in. 
d 3 = 0.25 in. 

/ = 16.7 in. 
t s = 0.06 in. 

t c = 0.0015 in. 


which gives a honeycomb core density of P 3 / P, = 0.019976 (face sheets excluded; eq. (1)) and a sandwich panel 
density of P^J/ Pc = 0.043497 (core plus face sheets; eq. (2)). For identical effective densities (that is, 
Pi = P 2 - P 3 - Pa = 0.019976P C or p[]] = p{^j = p{^} = p[^J = 0.043497 P c ), the dimensions of the four 
types of honeycomb cells considered are listed in table 1 , where for the corrugated cell, the range of corrugation 
depths h c /t c = 0-6 covers h c /t c = 0, 0.50, 1.00, 1.50, 1.75, 2.00, 2.25, 3.00, 3.50, 4.00, 4.50, 5.00, 5.50, 6.00. 
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Table 1. Dimensions of different honeycomb cells with identical effective core densities 
(Pi = p 2 = p 3 = p 4 = 0.019976P C ). 


Cell type 

h v 

b\ 

0(5 


*c ’ 

h c /t c 

P, 

(0 

in. 

in. 

in. 

in. 

in. 


in. 

Right hexagonal 
(/=1) 

0.1156 

0.1156 

0.2312 

0.0015 




Square 

a =2) 

0.1502 

0 

0.2124 

0.0015 




TPS (flat) 
O' = 3) 

0.1326 

0.0625 

0.2500 

0.0015 




TPS (corrugated) 
0 = 4) 

0.1326 

0.0625 

0.2500 


eq (12) 

0-6 

0.0265 


FINITE-ELEMENT ANALYSIS 

This section describes the heat-transfer analysis of unit honeycomb cells of different geometry, thermal 
bending analysis of a honeycomb TPS panel, and core-depth-wise buckling analysis of honeycomb cell walls of 
different geometry. The TPS panel is assumed to be fabricated with superalloy Inconel 617 (Inco Alloys 
International, Inc., Huntington, West Virginia). The Structural Performance and Resizing (SPAR) finite-element 
computer program (ref. 10) was used in the three kinds of analyses. 

Heat-Transfer Analysis 

The heat-shielding performance of the honeycomb TPS panel can be studied by simply considering the heat 
transfer in a single honeycomb cell. Figure 9 shows the finite-element thermal models generated for the four types 
of honeycomb cells for the heat-transfer analysis. The upper and lower face sheets were modeled with K31 
(ref. 10) triangular conduction elements with thickness t , and the honeycomb cell walls were modeled with 
quadrilateral K41 (ref. 10) conduction elements. Because each cell wall panel is shared by two adjacent cells, t / 2 
was used for the heat- transfer analysis as the thickness of K41 elements lying in the free-cell wall regions, and t 
for the thickness of K41 elements lying in the interbonded double-cell wall regions. 

For internal radiation exchanges, three-nodes R31 (ref. 10) radiation elements and four-nodes R41 (ref. 10) 
radiation elements were attached respectively to the inner faces of K31 and K41 elements. The internal radiation 
view factors were computed using a VIEW program (ref. 10) incorporated into the SPAR program. Because the 
honeycomb cell is an enclosed structure, the “ENCLOSURE” command was used to correct the calculated 
radiation view factors so that the sum of the view factors from each element is equal to unity. Namely, if 
Fjj ( i,j =1,2, 3,...) is the radiation view factor defined as the fraction of radiant heat leaving radiation element i 
and incident on radiation element j, then the “ENCLOSURE” command will impose the following condition for 
the final values of F ;i \ 




( 13 ) 


5=1 

j 


where the summation is taken over j for each given i. 

Further, the command “UNOBSTRUCT”; “ALL” 
(or radiation elements) are nonobstructive. Reference 
Inconel 617 for input to the thermal models. Table 2 
honeycomb cells. 


was used to define that all the internal radiation surfaces 
11 provided the thermal property data of the superalloy 
shows the sizes of thermal models for the four types of 


Table 2. Sizes of finite-element thermal models generated for 
different honeycomb cells. 


Cell type (i) 

JLOC 

K31 

K41 

R31 

R41 

Right hexagonal (/ = 1) 

525 

48 

480 

48 

480 

Square (i = 2) 

441 

40 

400 

40 

400 

TPS (flat) (/ = 3) 

525 

48 

480 

48 

480 

TPS (corrugated) (i = 4) 

525 

48 

480 

48 

480 


For the heat input to the upper face sheet of the honeycomb cell, a heating rate of 16 °F/sec was used. This 
heating rate is comparable to the steepest stage of surface heating rate at a certain body point of a new type of STS 
during its mission. Because the surface heating is a temperature input instead of a heat-flux input, external 
radiation exchanges with space could be neglected. Furthermore, because of the tiny volume inside each 
honeycomb cell, both convection and conduction heat transfer of the interior air mass were neglected. This section 
studies the effect of honeycomb cell geometry on the heat-shielding performance of the TPS panel. 

Thermal Bending Analysis 

Figure 10 shows a quarter-panel finite-element model generated for the honeycomb TPS panel for thermal 
bending analysis. Because of symmetry with respect to the in-plane x- and y-axcs, only one-quarter of the TPS 
panel was modeled. The commands “SYMMETRY PLANE = 1” and “SYMMETRY PLANE = 2” (ref. 10) were 
used to create the mirror images of the quarter-panel region, thereby simulating the whole panel. 

Each face sheet is modeled with one layer of quadrilaterally combined membrane and bending E43 (ref. 10) 
elements, and the honeycomb core with a single layer of hexahedron (or brick) S81 (ref. 10) elements connecting 
the upper and lower face sheet E43 elements. 

In the thermal bending analysis, the extensional and bending stiffness of the sandwich panel is assumed to be 
carried by the two face sheets, and the transverse shear stiffness by the honeycomb core. The TPS panel will be 
supported under the following different boundary conditions for comparing different degrees of thermal 
deformations and thermal stress fields induced in the TPS panel: 
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1. Fixed simply supported comers: four corners of the middle plane are simply supported and fixed. 

2. Free simply supported corners: four corners of the middle plane are simply supported, but are allowed to 
move freely in the x and y directions. 

3. Fixed clamped corners: clamped four corners are fixed. 

4. Free clamped corners: clamped four corners are allowed to move freely in the x and y directions. 

5. Fixed simply supported edges: four edges of the middle plane are simply supported and are fixed. 

6. Free simply supported edges: four edges of the middle plane are simply supported and are allowed to move 
freely in the x and y directions. 

7. Fixed clamped edges: clamped four edges are fixed. 

8. Free clamped edges: clamped four edges are allowed to move freely in the x and y directions. 

This method of simply supporting a sandwich panel is slightly different from simply supporting a solid flat plate 
(ref. 12). To simulate the fixed or free simply supported condition for the sandwich panel, pin-ended rigid rods were 
attached vertically to the panel edges (or comers) to connect the two face sheets (figs. ll(a)-(b)). Because the 
honeycomb core does not carry extensional or bending stiffness, edge support points cannot be attached to the 
honeycomb core. The midpoint of each rigid rod was pin-joined to a point (fixed or movable in the x-y plane) lying in 
the hypothetical panel middle plane (ref. 12). Each pin-ended rigid rod was modeled with two identical E22 (ref. 10) 
elements (beam elements for which the intrinsic stiffness matrixes are given). To simulate the rigidity of the rods, the 
extensional and the transverse shear stiffness of the E22 elements were made very large. The pin-joint condition at the 
face sheet edges was simulated by assigning zero values to the rotational spring constants in the stiffness matrix for the 
E22 elements. To simulate the pin-joint condition at the hypothetical middle plane, the three rotational constraints were 
eliminated. One node of each E22 element was connected to the associated node of E43 element, and the other node 
was connected to the hypothetical middle-plane point. For the fixed or free simply supported corners, only one pair of 
E22 elements was attached to one corner of the quarter-panel model. 

For the fixed or free clamped edges (figs. 1 1 (c)— (d)), the panel edges (or corners) were constrained to have 
zero transverse rotations. The sizes of the quarter-panel models for the honeycomb TPS are given below. 


Table 3. Sizes of quarter-panel finite-element 
models for honeycomb TPS panel. 


Support condition 

JLOC 

E22 

E43 

S81 

Simply supported, 

1299 

98 

1152 

576 

clamped edges 





Simply supported. 

1251 

2 

1152 

576 

clamped corners 






The upper face sheet will be heated uniformly to temperature, T u , of 800 °F and the lower sheet to temperature, 
Tj, of 450 °F. These temperature levels are similar to the temperatures across the TPS depth giving the maximum 
temperature gradient at a certain body point during a new STS mission. The material properties of the Inconel 617 
were obtained from reference 11, and the effective elastic constants of Inconel 617 honeycomb core (with TPS cell 
geometry) were calculated using the formulae developed in reference 13. This section examines the panel 
deflections and the levels of the induced thermal stresses in the TPS panel under different support conditions. 
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Cell Wall Buckling Analysis 


Compressive buckling behavior of a honeycomb panel subjected to surface aerodynamic pressure q (fig. 1) can 
be studied by simply considering the weakest cell wall panel of each type of honeycomb cell. Figure 12 shows the 
locations of the weakest cell wall panels analyzed, and figure 13 shows the corresponding finite-element models. 
For simplicity, the cell wall panel boundaries are simply supported. For the corrugated cell wall panel, two cases 
are considered: varying thickness ~t c to maintain constant density (P 4 = P 3 ); and constant panel thickness t (to 
determine how p 3 increases with increasing h ( . ). 

Let <x be the compressive stress in the honeycomb cell walls induced by the surface pressure q. At the 
buckling of type i cell walls (that is, 0 , = ( 0 J '.]. ), the critical surface pressure q = q^J can be calculated from 
the following: 


AO 

(0 = _jv_, M) 
“cr (iV U z'cr 

K 


P,^zt 


(14) 


(4) 

where, for the corrugated cell panel, the panel load (AM is associated with the unit length along the 

~ cv M 

coiTugation length, and for the constant density case (P 4 = P 3 ) , the panel thickness t is replaced with t c . 

Applying the expressions for the effective densities P,/ P c (eqs. (4)-(7)), the buckling pressure q^j 
associated with type i cell can be written as follows: 


Type 1: 


Type 2: 


Type 3: 


Type 4: 


U> = I6J3. \( 1 ) 

' 9rf, ' z’er 


P = U \ N )< 2 > 

' z’er 


(3) = 

q ° r 15 d 3 


(^ + 3 )(^ 3 ; 


(15) 


(16) 


(17) 


(4) = %J2_ 
lcr 15d 4 


/ 

72 + 3 

\ 



(ap) (4) 


(18) 
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The panel buckling loads (/V,)^, (N , (N S)^ for the flat cell wall panels of honeycomb cells types 1, 2, 

and 3 can be easily calculated from the classical buckling equation for a solid flat rectangular plate (length a, width 
b t ) (ref. 14): 


(N ) (0 
z'cr 


kD / 


2 2 \ 


, n a 

m H 

2 


m 


b • 


i = 1, 2, 3 


z / 


(19) 


However, finite-element buckling solutions for the flat cell wall panels of honeycomb cells types 1, 2, and 3 were 
also obtained and compared with a Timoshenko solution (ref. 14) (eq. (19)). The accuracy of those solutions give 
confidence in the finite-element buckling solutions obtained for the corrugated panel using similar element density. 

For simply supported slender rectangular panels like those with cell wall types 1, 2, and 3, the panels will 
buckle into several half-waves (m * 1 ) in the loading direction, and only one half-wave ( n = 1 ) in the transverse 
direction, (ref. 14). Taking n = 1 in equation (19) and minimizing (N z )^j by differentiating equation (19) with 
respect to m, and then setting the resulting equation to zero, leads to 


m 


b ■ 


(20) 


which gives the correct value of m if a/b- is an integer. If a/b- is not an integer, the correct value of m should be 
an integer closest to the value of a/b . . 

Alternatively, the value of m for minimum buckling load of a rectangular panel with an aspect ratio of a/b-, 
can be found graphically from figure 14, the plots of a family of buckling curves associated with different m for 
rectangular plates (ref. 14). Namely, first locate the value of a/b- on the abscissa, and draw a vertical line from 
point a/ b j . This vertical line will intersect the curve segment of the lowest buckling load (solid curves, fig. 14), 
and thereby find the value of m associated with that particular buckling curve segment. Table 4 lists the aspect 
ratios a/b i and the associated m for minimum buckling loads determined for the flat wall panels of types 1, 2, and 
3 cells. 


Table 4. Aspect ratios and number of buckled half-waves 
for honeycomb cell wall panels. Pi = P 2 - P3 • 


Cell type (, i ) 

a/bj 

m 

n 

Right hexagonal (i - 1) 

4.22 

4 

1 

Square (/' = 2) 

3.25 

3 

1 

TPS (flat) (, i = 3) 

3.68 

4 

1 
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RESULTS 


The effect of honeycomb cell geometry on the heat-shielding performance and cell wall buckling 
characteristics of a superalloy TPS has been investigated. The results of heat-transfer, thermal bending, and cell 
wall buckling analyses of superalloy honeycomb TPS panels are presented in the following sections. 

Heat Transfer 

Figure 15 shows the linear time history of the temperature input to the outer face sheet (point A), and the time 
histories of the temperatures at the inner face sheet (point B) for different cell geometry. The difference between 
the outer (T u ; point A) and the inner ( T / ; point B) face sheet temperatures. AT\ (T - T,) , is the measure of the 
heat-shielding performance of the TPS. Namely, the larger the value of AT , the better the heat- shielding 
performance. The AT for all cell geometries reaches maximum at approximately 50 sec, then decreases only 
slightly with the increasing time, t. The temperatures at point B for all cell geometries are nearly the same, 
indicating that the TPS heat-shielding performance is relatively insensitive to the shape change of the honeycomb 
cell (under the same effective density). The right hexagonal cell has the lowest heat-shielding performance 
(the lowest AT), and the square cell has the highest heat-shielding performance (the highest AT). The temperature 
curves at point B for the two types of TPS cells lie between those for the right hexagonal and square cells. For 
conduction only (without internal radiation), all the temperature curves at point B for the four types of cell 
geometry are almost identical, and are hardly distinguishable. The effect of internal radiation turned out to be much 
smaller than that of conduction for the present TPS core geometry ( a = 0.488 in.; d 3 = 0.250 in.). For example, 
at t = 100 sec (fig. 15), the temperature at point B for the TPS cell (flat wall) consists of 13-percent radiation 
component and 87-percent conduction component. Table 5 shows typical data of figure 15 at t = 50 sec and 
t = 100 sec where the “heat-shielding index” is defined as the ratio of AT of any cell geometry to the AT of the 
TPS cell with flat cell walls (/ = 3) which was used as the reference cell geometry. 


Table 5. Fleat-shielding performance of superalloy TPS panels with different honeycomb-core 


cell geometry. P\ - P 2 = 

P 3 = Pa 

= 0.019976P C . 





t = 50 sec 


t = 100 sec 

Cell type ( i ) 

AT, °F 

Fleat-shielding index 

AT, °F 

Fleat-shielding index 

Right hexagonal ( i = 1) 

432 

0.99 

355 

0.95 

Square (i = 2) 

443 

1.02 

409 

1.09 

TPS (h c = 0) (i = 3) 

436 

1.00 

375 

1.00 

TPS (h c /t c = 6) (i = 4) 

440 

1.01 

391 

1.04 

Reference cell geometry. 
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Figure 16 shows the effects of modified core depth a' and cell size d 3 on the heat-shielding performance of 
TPS with flat cell walls. For the core depth a' = a , changing the cell size from r/ 3 = 0.25 in. to d 3 = 0.5 in. 
(reducing effective density) increased AT slightly during the initial heating stage until t = 65 sec, but reduced 
AT after t = 65 sec. For the deeper core depths (o' = 1.5o and o' = 2 o ), the honeycomb core with the smaller 
cell size (d-, = 0.25 in.) has better heat-shielding performance (larger AT) than those with the larger cell size 
( d-, = 0.5 in.), which have lower effective density. The comparatively lower values of AT for the larger cell size 
d 3 = 0.5 in. is caused by the increased radiation effect from the upper face sheet to the lower face sheet. Table 6 
shows typical data from figure 16 at t = 50 sec and t = 100 sec, where the “weight index” is defined as the ratio of 
the weight of an / x / TPS panel with any cell size r/ 3 and depth o' to the weight of an / x / TPS panel with cell 
size r/ 3 = 0.25 in. and depth o' = o = 0.488 in. calculated using equation (3). 


Table 6. Fleat shielding performance of TPS honeycomb-core with different cell 
sizes and core depths; flat cell walls (/ = 3). 


Cell size 
(d 3 ), in. 

a ' , 
in. 


t = 50 sec 


t = 100 sec 

Weight 

index 

AT 

Fleat- shielding 
index 

AT 

Fleat-shielding 

index 

0.250 

0.488 (In) 

436 

1.00* 

375 

1.00* 

1.00* 

0.250 

0.732 (1.5a) 

592 

1.36 

633 

1.69 

1.25 

0.250 

0.976 (2a) 

682 

1.56 

886 

2.36 

1.50 

0.375 

0.488 (la) 

455 

1.04 

321 

0.86 

0.83 

0.375 

0.732 (1.5a) 

592 

1.36 

531 

1.42 

1.00 

0.375 

0.976 (2a) 

674 

1.55 

754 

2.01 

1.17 

0.500 

0.488 (la) 

466 

1.07 

284 

0.76 

0.75 

0.500 

0.732 (1.5a) 

587 

1.35 

454 

1.21 

0.87 

0.500 

0.976 (2a) 

664 

1.52 

642 

1.71 

1.00 


Reference cell geometry. 
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Thermal Bending 


Figures 17 and 18 respectively show typical deformed shapes of honeycomb TPS panel subjected to thermal 
bending (T = 800 °F; T l = 450 °F) for cases with simply supported fixed corners and simply supported fixed 
edges. Table 7 shows a summary of the deflections at the sandwich panel edge midpoint, w g , and center, w max , 
and stresses { a , cr } induced at the centers of the upper and lower face sheets. 

x y 

Note that only for the “free simply supported corners” case (marked with an asterisk in table 7) are the face 
sheets stress-free; however, the TPS panel bulged out 0.4144 in. To eliminate TPS outward bulging, the panel 
support condition must be either “fixed clamped edges” or “free clamped edges” (marked with a dagger in table 7). 
These two support conditions, however, will certainly induce higher stresses in the face sheets. The simply 
supported and clamped fixed edge cases induced extremely high face sheet stresses that could cause the entire TPS 
panel to buckle (ref. 15). 


Table 7. Panel deflections and thermal stresses induced in the face sheets of superalloy honeycomb TPS under 
different support conditions; T u = 800 °F ; Tj = 450 °F ; TPS cell (flat walls, i = 3). 


Boundary 

condition 

Panel 

center 

Edge 

midpoint 

Upper 

center 

Lower 

center 

w max ’ i n - 

w e > in- 

a x , lb/in 2 

o y , lb/in 2 

o x , lb/in 2 

o y , lb/in 2 

Fixed simply supported 

0.4144 

0.2061 

-25,060 

-25,060 

-25,060 

-25,060 

corners 







Free simply supported corners* 

0.4144 

0.2061 

0 

0 

0 

0 

Fixed clamped corners 

0.2926 

0.1085 

-33,820 

-33,870 

-16,300 

-16,240 

Free clamped corners 

0.2926 

0.1085 

-8,760 

-8,820 

8,760 

8,820 

Fixed simply supported edges 

0.1689 

0 

-201,420 

-201,360 

-166,340 

-166,410 

Free simply supported edges 

0.1689 

0 

-17,540 

-17,480 

17,540 

17,480 

Fixed clamped edges ' 

0 

0 

-241,390 

-241,390 

-126,370 

-126,370 

Free clamped edges ’ 

0 

0 

-57,510 

-57,510 

57,510 

57,510 


^Support condition for stress-free face sheets. 
^Support conditions for zero TPS panel deflection. 
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Cell Wall Buckling 


Figure 19 shows buckling shapes of the four types of honeycomb cell wall panels under axial compression 
caused by air load. As predicted in table 4, the right hexagonal cell wall and the TPS cell flat wall buckled into four 
half-waves in the cell axial direction (figs. 19(a) and 19(c)). The square cell wall panel buckled into three 
half-waves (fig. 19(b)) in the cell axial direction; and the TPS cell corrugated wall ( h/t = 6 ) buckled into only 
one half-wave in the cell axial direction because of corrugation (fig. 19(d)). Table 8 shows a summary of the 
numerical results of honeycomb cell wall buckling analysis, and that the buckling solutions obtained from the 
SPAR finite-element method and Timoshenko classical elasticity theory (ref. 14) are extremely close, giving 
confidence in the accuracy of the SPAR buckling solution obtained for the corrugated cell wall case (/ = 4) . 


Table 8. Buckling pressures of superalloy honeycomb TPS with different cell geometry. 


Cell type (i) 

q l c l , lb/in 2 

Timoshenko SPAR 

Solution difference, 
percent 

q«W 3) 

^ cr ^ cr 

Right hexagonal ( / = 1 ) 

214.7460 

214.8327 

0.0404 

1.31 

Square (i = 2) 

127.7003 

127.7508 

0.0395 

0.78 

TPS (h c = 0) (/ = 3) 

163.9307 

164.1705 

0.1463 

1.00* 

TPS (h/t c = 6) (i = 4) 

No theory 

682.0951 

N/A 

4.16 


Reference cell geometry. 


As shown earlier, the heat-shielding performances of different types of honeycomb cells are quite similar. The 
air pressure buckling strengths, however, are quite different (table 8). The cell axial buckling strengths of the right 
hexagonal and square cells respectively are 1.31 and 0.78 times that of the TPS cell with flat walls. The TPS cell 
with coiTugated walls ( h c /t = 6) has a buckling strength 4.16 times that of the TPS cell with flat walls (h ( . = 0) . 

The maximum dynamic pressure a Space Shuttle orbiter has experienced during reentry is approximately 
1.84 lb/in (265 lb/ft") at Mach 0.6. Considering that the new STS will be subjected to similar dynamic pressure, 
these four types of honeycomb cells certainly have excess buckling strengths to carry the surface air load. 

Figure 20 shows plots of buckling pressure q'^ and core density p,-/ p c of type i cells as functions of 
honeycomb cell size d . . This figure provides vital design curves for the honeycomb cells. Figure 20 is used as 
follows: for a given fixed core density p,-/ p c = 0.019976 , a horizontal line ABC is drawn to intersect the core 
density curves p 3 , pj, and p 2 for cell types 3, 1, and 2. Then vertical lines are drawn from points A, B, and C. The 
intersection of those vertical lines with the abscissa then gives the cell sizes d 3 = 0.25 in., d^ = 0.2312 in., and 
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= 0.2124 in. The intersections of the vertical lines drawn from points A. B. and C with the q^ r curves give the 

values of buckling pressure q^_ associated with each cell type; namely, q'^ = 163.93 lb/in 2 , 
(11 i PI i 

q cr = 214.75 lb/in", and q cr = 127.70 lb/in". Note that the right hexagonal cell has the highest buckling 
strength, and the square cell has the lowest buckling strength. The value of buckling pressure of the TPS cell with 

(3) 

flat walls ( q cr ) lies between those of the right hexagonal cell and square cell (table 8). 

(4) (3) 

Figure 21 shows the plots of buckling pressure q cr /q cr , core density P4/P3, and sandwich density 
(4) (3) 

p sw / p sw of TPS cells with corrugated walls as functions of the corrugation depth h c /t c . Data for plotting the 
curves associated with this cell are normalized by the corresponding values associated with the TPS cell with flat 
walls (type 3 cell). Figure 21 shows two cases of buckling curves. The first case has a constant sandwich density 
(that is, p'^J = p[ 3 J ), but t is a function of corrugation depth hJt n . The second case has a constant corrugation 

a W IV c C L 

(4) (3) 

wall thickness (that is, t c = t (; ), but the core density P4/P3 and the sandwich density p sw / p sw can change with 
the corrugation depth h c /t c . 

At the lower value of h £ / 1 ( . ( h /t < 0.5) , the corrugated panel behaved like the flat wall panel (table 4) and 
buckled into four half-waves (m = 4, n = 1) in the cell axis direction. The corrugated panel then buckled into 
(m = 3, n = 1) mode in the region 0.5 < (li r /t ( .) < 2, and then into (m = 2, n = 1) mode in the region 
2 < h c /t c < 5.5 , and finally into (m = 1, n = 1) mode when the coiTugation depth /q /? ( . was increased beyond 
h/t c = 5.5 . 

Buckling pressure q^/q^J of the second case (t ( . = constant) is slightly higher (by 1.83 percent at 
h c /t c -6) than that of the first case because in the first case, the cell wall thickness ~t c decreased with the 
increasing corrugation depth h c /t c and thus lowered the buckling pressure. By keeping the wall thickness 
constant, ~t c = t , independent of the corrugation depth change (case 2), the honeycomb core density p 4 / p 3 and 
sandwich panel density p^V p{ 3 j increased by merely 3.66 and 1.64 percent, respectively, at h /t r = 6. These 
amounts of density sacrifices could be considered insignificant. 

CONCLUSIONS 

Fleat-transfer, thermal bending, and mechanical buckling analyses were performed on a superalloy thermal 
protection system (TPS) for future hypersonic flight vehicles. Effect of honeycomb cell geometry on the 
heat-shielding performance and honeycomb cell wall buckling characteristics, and the effect of boundary 
conditions on the TPS thermal bending, were investigated. 

The results are summarized as follows: 

1. The heat- shielding performance of a honeycomb TPS is insensitive to the shape of the honeycomb cell 
under the same effective core density, but improves with the core depth. 

2. The TPS outward bulging effect can be eliminated if the panel support condition is either “fixed clamped 
edges” or “free clamped edges.” The latter case is preferred because of much lower stresses. 

3. To reduce the TPS face sheet stresses to zero, the boundary support condition must be “free simply 
supported corners.” 
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4. The cell axial buckling strength is somewhat sensitive to the shape of honeycomb cell, but is quite sensitive 
to the depth of cell wall corrugation. The right hexagonal cell and square cell have axial buckling strengths 
1.31 and 0.78 times, respectively, that of the TPS cell with flat walls; and the TPS cell with corrugated 
walls with corrugation depth h /t = 6 has buckling strength 4.16 times that of the TPS cell with flat 
walls. 

5. Through corrugation of the honeycomb cell walls, the buckling pressure of the TPS panel could be 
enhanced greatly. 


Dryden Flight Research Center 
National Aeronautics and Space Administration 
Edwards, California 
March 1, 2001 
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FIGURES 



Figure 1. Honeycomb sandwich thermal protection system (TPS) subjected to heating T(t) and aerodynamic 
pressure q over entire upper surface. 
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Figure 2. Formation of honeycomb structure through pulling apart of joined thin strip foil assembly. 
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Figure 7. Relationship between TPS cell wall corrugation angle, 6, and corrugation depth h /t . 



Figure 8. Decreasing of TPS cell corrugated wall thickness t c /t with the increasing corrugation depth V'c 
under constant core density; p 4 = p 3 . 
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(a) Right hexagonal cell. 


(b) Square cell. 



(c) TPS cell (flat walls). 



(d) TPS cell (corrugated walls). 


Figure 9. Finite-element models for different unit honeycomb cells for heat transfer analysis. 
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(a) Fixed simply supported edge or corner. (b) Free simply supported edge or corner. 
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(c) Fixed clamped edge or corner. (d) Free clamped edge or comer. 


Figure 1 1 . Different support conditions at edges or corners of honeycomb sandwich TPS panel. 
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(a) Right hexagonal cell wall panel. 





Figure 12. Different honeycomb unit cell wall panels for compressive buckling analysis. 
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(a) Right hexagonal cell wall panel. 


(b) Square cell wall panel 



(d) TPS cell corrugated wall panel. 


Figure 13. Finite-element models for different honeycomb cell wall panels for compressive buckling analysis 
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Figure 16. Effects of cell sizes and core depths on the heat-shielding performance of superalloy honeycomb TPS 
panel; TPS cell with flat cell walls. 
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Figure 17. Deformed shape of honeycomb sandwich TPS panel under differential surface heating; T 
T , = 450 °F; simply supported fixed corners. 


Figure 18. Deformed shape of honeycomb sandwich TPS panel under differential surface heating; T 
T, = 450 °F; simply supported fixed edges. 
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(a) Right hexagonal cell wall panel. 
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(b) Square cell wall panel. 
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(c) TPS cell flat wall panel. 


(d) TPS cell corrugated wall panel. 


Figure 19. Buckled shapes of different honeycomb cell wall panels; simply supported. 
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Figure 21. Plots of buckling pressure / q^J and core density p 4 /p 3 of corrugated TPS cell as functions of 

corrugation depth h /t . 
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APPENDIX 


HONEYCOMB CELL ANALYSIS 


The geometrical analyses of the four types of honeycomb cells are performed below. In each honeycomb cell 
are two interbonded double cell walls of thickness 2 1 (interbonded region diminished theoretically for the 
square cell), and four free cell walls with single thickness t ..In the calculations of the wall cross-sectional area 
A u .' of the unit cell, only one-half of cell wall thickness is used because each cell wall panel is mutually shared 
by two adjacent cells. Namely, for the free cell walls, thickness t /2 is used; and for the bonded double walls, 
thickness 2t c /2 = t is used. 


Right Hexagonal Cell: 


From figure 3(a), the length b j of one side of the right hexagon is related to the cell size d x through 


b \ 



(A-l) 


The cross-sectional area A 
calculated as 


(1) 

c 


of the unit right hexagonal cell (consisting of six right triangles) can be 



= 6 


~2 

\ / 


&b l 

2 1 




373,2 
2 1 


(A-2) 


The cross-sectional area of the unit hexagonal cell wall, which consists of two bonded sides and four free 
sides, can be calculated as follows: 


A 


(1) 

w 



+ 4 b x 



4b l 7 


(A-3) 


From equation (A-2) and (A-3), the effective density P\ of the right hexagonal cell can be expressed as 


Pi _ _ 873 t c _ 1673 t c _d x 

P c A 0) 9 b x 9d x ’ 1 2 


(A-4) 


Square Cell: 

From figure 3(b), the length b 1 of one side of the square cell is related to the cell size as 
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b = 


J2d~, 


(A-5) 


( 2 ) 

The cross-sectional area A of the unit square cell is given by 


A 2> - A 


(A-6) 


( 2 ) 

The cross-sectional area A of the square cell wall (consisting of four free sides) is expressed as 


- 4b 


(t c \ 

A I = 2 b*t 


2 c 


(A-7) 


In light of equations (A-6) and (A-7), the effective density P 2 of the square cell is given by 


^ = = ; b, = ^ 

Pc b 2 d 2 2 2 


(A-8) 


TPS Cell (Flat Wall): 


From figure 3(c), the length b~, and b' are respectively related to the cell size d-> as 


b 3 = 


_ 2>j2, . , , _ ^3 _ J2. 


= -r = ^b. 


(A-9) 


( 3 ) 

The cross-sectional area A , of the unit TPS cell (consisting of two identical trapezoids) can be calculated as 


A 31 = 2 


Ifdj + b') 

> J2. 


5,2 


= ^ b 


(A- 10) 


(3) 

in which equation (A-9) was applied. And, the cross-sectional area A of the unit TPS cell wall (consisting of 
four free sides and two bonded sides) can be calculated as 


A 3) = 


it s \ 

V 2 / 


+ 2 b' 


(2t c \ 


\ 2 ) 


= 2 (b, + b')t r 


(A- 11) 


which becomes, in light of equation (A-9), 


A 3, = !(3 + V2)Vc 


(A- 12) 
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From equations (A-10) and (A-12), the effective density p 3 of the TPS cell (flat wall) can be expressed as 


^(3) 

fh K 2 


Pc A 


(3) 5 


= if 3 + 72)? = ^(3 + 74? ; 6, = 


63 15 


/J3 ’ ' 3 


_ 3V2 


(A- 13) 


TPS Cell (Corrugated Wall): 


From figure 3(d), the length b 4 and b' arc respectively related to the cell size d, as 


_ 3j2 


_ 72, 


b, = ^L± d ; b' = -= = ^±b, 
4 § 4 4 3 4 


(A- 14) 


(4) 

Considering equation (A-9), the cross-sectional area A ( , of the unit TPS cell (consisting of two identical 
trapezoids) can be calculated as 


4 4) = 2 


i{d A +b\ b * 


> Jl 


5,2 
= ^4 


(A- 15) 


( 4 ) 

And, the cross-sectional ai ca A w of the unit TPS cell wall (consisting of four free sides and two double sides) 
can be calculated as follows: 


C 41 = 46 


1 

2 

\ J 


+ 2 b' 




2 

V / 


- 2 [b + b'Yt 


(A- 16) 


which becomes, in light of equation (A- 14), 


A (4) = 2 


i + 4 \i c 


(A-17) 


From equations (A-15) and (A-17), the effective density p 4 of the TPS cell (corrugated wall) can be expressed 


as 


P_4 

Pc 




(A- 18) 


where the corrugated side length b is calculated in the following. 


The TPS cell with corrugated wall (fig. 3(d)) has identical size and cell density as the TPS cell with flat wall 
(fig. 3(c)). Therefore, to maintain constant effective density P4 = p 3 , the corrugated wall thickness ~t c must be 


36 



modified whenever the corrugation depth h is changed. From figure 4, the corrugation true length b consists of 
five identical corrugation legs of length c. Namely, 

b = 5c (A- 19) 

The length c of each corrugation leg is related to the corrugation half-pitch p and the corrugation depth h ( , 
through 


c = Jp 2 + K 


(A-20) 


where the half-pitch length p is related to b, (end-to-end distance of con'ugation side (fig. 4)) and h through 


5 P = 


(A-21) 


Combining the above three equations results in 


b_ 
b , , 




1 + 24— = 




, 256 h c r 3 V2 , 

l + ~j2 ' * 4 = — ^4 
C* 4 


(A-22) 


Substituting equation (A-22) into equation (A- 18) yields 
4(4) 


P4 

Pc 


w 

(4) 

c 


/ 

3 

h] r ) 

1 + 24 — + 72 

tc _ 8 72 

t 

3 

ll + 256^£ + 72 

1 \ 

K J 

b 4 15 

i \ 

3 A J 


— ■ b - 

A ’ - 


_ 372 


(A-23) 


For the constant effective density (P 4 = P 3 ) (that is, = A^), the con'ugated cell wall thickness t c is 
related to t , through 


2 {b + b')tc = 2(b 3 + b'\t c 


(A-24) 


which is rewritten as 


tc = b 3 + b ' 

7 b + b' 


(A-25) 


which, through the application of bn = b 4 and equations (A- 14) and (A-22), leads to 


tc 

t„ 


3 + 72 


3 + 72 


ft c r 

1 + 24— + 72 3 

2 


; b A = 


372 


K r 
1 + 256— + 72 

d 4 


(A-26) 
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Substituting equation (A-26) into equation (A-23) yields 


p 4 _ A ( J _ 2(3 + J~2)t c _ 872(3 + J2)t c _ 2>/ J~ 2 
P c /*) 5 b A \5d A ’ 4 8 

which is si mi lar to equation (A-13) (the flat cell wall case, i = 3 ). 

NOTE 


(A-27) 


This report stemmed from the investigations of heat shielding and thermoelastic performance of superalloy 
thermal protection system (TPS) to be used for protecting the X-33 research vehicle from overheating during its 
mission. The report contains no classified information, and serves as a stimulant for university professors and their 
graduate students to extend Ko's basic research work. 

William L. Ko 
April 27, 2000 


Material Properties 

Table A-l. Material properties of Inconel 617 used for superalloy TPS panel. 


T, °F 

450 

625 

800 

1800 

2000 

E, lb/in 2 

28.25 x 10 6 

27.50 x 10 6 

26.50 x 10 6 

18.20 x 10 6 

16.0 x 10 6 

G, lb/in 2 

1 1.00 x 10 6 

10.60 x 10 6 

10.20 x 10 6 

6.60 x 10 6 

5.71 x 10 6 

V 

0.284 

0.297 

0.300 

0.390 

0.400 

a, in/in- °F 

7.12 x 10" 6 

7.50 x 10“ 6 

7.65 x 10" 6 

9.10 x 10“ 6 

9.25 x 10" 6 

p c , lb/in 3 

0.3015 

0.3015 

0.3015 

0.3015 

0.3015 
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